A Parametrization of Bipartite Systems Based on S77(4) Euler Angles 



Todd Tilma 

The Ilya Prigogine Center for Studies in Statistical Mechanics and Complex Systems 

Physics Department 
The University of Texas at Austin 
Austin, Texas 78712-108^ 

Mark Byrd 
Maxwell Duiorkm Laboratory 
33 Oxford Street 
Harvard University 
Cambridge, Massachusetts 0213^ 

E.C.G. Sudarshan 
Center for Particle Physics 
Physics Department 
The University of Texas at Austin 
Austin, Texas 78712-108^ 
(Dated: February 3, 2008) 

In this paper we give an explicit parametrization for all two qubit density matrices. This is 
important for calculations involving entanglement and many other types of quantum information 
processing. To accomplish this we present a generalized Euler angle parametrization for SU(4) 
and all possible two qubit density matrices. The important group-theoretical properties of such 
a description are then manifest. We thus obtain the correct Haar (Hurwitz) measure and volume 
element for SU(4) which follows from this parametrization. In addition, we study the role of this 
parametrization in the Peres-Horodecki criteria for separability and its corresponding usefulness in 
calculating entangled two qubit states as represented through the parametrization. 



* Email: ;ilma@Dhvsics.utcxas.edv. 

t Email: :>rd@hrl. harvard. cdtj 

t Email: 3udarshan@physics.utcxas.edu 



2 



I. INTRODUCTION 



In quantum mechanics the appropriate description of mixed states is by density matrices. For example, their 
compact notation makes them useful for describing entanglement and decoherence properties of multi-particle quantum 
systems. In particular, two two-state density matrices, also known as two qubit density matrices, are important for 
their role in explaining quantum teleportation, dense coding, computation theorems, and other issues pertinent to 
quantum information theory. 

Although the ideas behind extending classical computation and communication theories into the quantum realm 
have been around for some decades now, the first reference to calling any generic two-state system a qubit comes 
from Schumacher [0 in 1995. By calling a two-state system a qubit, he quantified the relationship between classical 
and quantum information theory: a qubit can behave like a classical bit, but because of the quantum properties 
of superposition and entanglement, it has a much larger information storage capacity. It is this capacity to invoke 
quantum effects to increase information storage and processing, that gives qubits such a central role in quantum 
information theory. 

Now, a qubit is just a state in a two-dimensional Hilbert space ||. If H ~ C 2 in the vector space, the unit vectors 

|^)=o|0) + 6|l) (1) 

with a and b being complex numbers satisfying 

M 2 + |6| 2 = l (2) 

define, up to a phase, the pure quantum states. In quantum information theory, the orthonormal basis {|0) , |1)} is 
used to represent the bit states (off) and 1 (on). As pointed out by Brown || the physical representation of these 
two bit states depends on the "hardware" being discussed; the basis states may be polarization states of light, atomic 
or electronic spin states, or the ground and first excited state of a quantum dot. 

If the qubit represents a mixed state, which is quite often the case, one should use a two-dimensional density matrix, 
which was introduced independently by Landau and von Neumann in the 1920's (see for example the discussion in 
) , for its representation. The formalism of density matrices allows one to exploit simple matrix algebra mechanisms 
to evaluate the expectation value of any physical quantity of the system. More recently, it has been pointed out by 
several people (see [|[ |[ |6[ [z|,and references within) that the density matrix representation of quantum states is also 
a very natural representation to use with regards to quantum information calculations. 

Following this lead we therefore express one qubit as 

p=^(h+cr-n), (3) 

i.e. as a general 2 by 2 hermitian matrix with unit trace and the positivity condition Tr[p] > implying n ■ n < 1 
or p 2 < p. Therefore, these density matrices are the disk D 3 , whose boundary <9L> 3 = S 2 = CP 1 represents the pure 
states (p 2 = p, or n ■ n = 1), and which can be thus characterized by the two angles < 9 < ir (the latitude) and 
< <\> < 2ir (the longitude) of the sphere S 2 . 

Now, two-state density matrices live in a 3 by 3 hermitian-matrix space, with CP 2 — SU(3)/U(2) as a subspace 
of pure states. Much is already known about these two- and three- state density matrices, especially when one uses, 
for example, Euler angle parametrizations (see |6) for more information). But what is not as well known is how the 
density matrices of larger dimensional Hilbert spaces, and thus of multiple qubits, looks under such a parametrization. 
This paper will make a great deal of progress in remedying this situation by giving an explicit parametrization of 
the density matrix of two qubits that is not redundant in the representation of the corresponding four-dimensional 
Hilbert space, and at the same time offers up the natural (Bures) volume measure on the set of all two qubit density 
matrices. This will be achieved by starting with a diagonal density matrix p,i, which represents our two qubit system 
in some particular basis, and then performing a unitary (U^ 1 = W), unimodular (Det[C/] = 1) transformation 

Pd -> U Pd U* (4) 

for some U £ SU(A), thus describing p in an arbitrary basis g, |[ ||, [To|. Wc should point out that progress in this 
direction has also been made in a recent publication by Schlienz and Mahler 111). 



II. EULER ANGLE PARAMETRIZATION FOR SU(4) 



We begin by giving the Euler angle parametrization for SU(4). Define U G SU(4). Using the Gell-Mann basis for 
the elements of the algebra (found in Appendix |X|) , the Euler angle parametrization is then given by 

„iA3»l iA2Q2 iA 3 a 3 iA 5 a 4 iA 3 a 5 c i\ioas iA 3 a 7 iA 2 Q8 JA 3 ag o iA 5 ai jA 3 an i= iA 2 ai2 iA 3 ai 3 ^iAsa^ iAi 5 ai 5 



3 



The derivation of this result is as follows. We begin by following the work of Biedenharn [12] and Hermann fl3] in order 
to generate a Cartan decomposition of SC/(4). First, we look at the 4 by 4, hermitian, traceless, Gell-Mann matrices 
Xi. This set is linearly independent and is the lowest dimensional faithful representation of the SU(A) Lie algebra. 
From these matrices we can then calculate their commutation relations, and by observation of the corresponding 
structure constants fijk we can see the relationship in the algebra that can help generate the Cartan decomposition 
of SU(A) (shown in detail in Appendix 

We now establish two subspaces of the SU (4) group manifold hereafter known as K and P. From these subspaces, 
there corresponds two subsets of the Lie algebra of SU{4), L(K) and L(P), such that for k\, fc 2 S L(K) and pi,P2 £ 
L(P), 

[k u k 2 ]£L(K), 
[pi,p 2 ] €L(K), 

[k l7 p 2 ]eL(P). (6) 

For 5/7(4), L(K) = {Ai, . . . , As, A15} and L(P) — {A9, . . . , A14}. Given that we can decompose the 5/7(4) algebra 
into a semi-direct sum |14{ 

L(SU(i))=L(K)®L(P), (7) 
we therefore have a decomposition of the group, 

U = K P. (8) 

From |l5| we know that L(K) contains the generators of the 5/7(3) subalgebra of 5/7(4), thus K will be the /7(3) 
subgroup obtained by exponentiating the subalgebra, {Ai, . . . , As} combined with A15 and thus can be written as (see 
H H for details) 

K — e i^3a e i\ 2 p e i\ 3 -t e i\ 5 6 e i\ 3 a e i\ 2 b & i\ 3 c e i\ & x e i^is<l> (g) 

Now, as for P, of the six elements in L(P) we chose the A2 analogue, A10, for 5/7(4) and write any element of P as 

P = K'e lX ^K" (10) 

where K' and K" are copies of K. 

Unfortunately, at this point in our derivation, we have a U with 28 elements, not the requisite 15 

U = KK'e iXl ^K". (11) 

But, if we recall that U is a product of operators in SU(4), we can "remove the redundancies," i.e. the first K' 
component as well as the three Cartan subalgebra elements of 5/7(4) in the original K component, to arrive at the 
following product ||, || 

JJ — e 'iA3Ql e iA2a2 e iA3Q3 e iA5Q4 e iA 3 a5 e iA2r; e iAio^ e iA3Q e iA2/3 e iA37 e iA5e e iA3a e iA26 e iA 3 c e iA s x e iAi50 ^2) 

By insisting that our parametrization must truthfully reproduce known vector and tensor transformations under 
SU{4), we can remove the last "redundancy," e lA2?? , and, after rewriting the parameters, generate equation (||) 

JJ _ giA 3 0!i gj'A 2 Q2 e iA 3 a 3 ^i\^a A ^1X305 giAi a6 giA 3 ct 7 ^i^ag ^iXaag ^i\ 5 a 10 e i\ 3 an ^i\ 2 ai2 e i^3a 13 ^iXsau ^i\ 15 a 15 (13) 

For our purposes it is enough to note that this parametrization is special unitary by construction and can be shown 
to cover the group by modifying the ranges that follow and substituting them into the whole group matrix, or into 
the parametrization of the characters fig] . 



III. DERIVATION OF THE HAAR MEASURE AND CALCULATION OF THE GROUP VOLUME FOR 

SU(4) 

Taking the Euler angle parametrization given by equation ([5]) we now wish to develop the differential volume 
element, also known as the Haar measure, for the group 5't/(4). We initially proceed by extending the method used 
in ||, ||] for the calculation of the Haar measure for SU(3); take a generic U £ SU(A) and find the matrix 



U- 1 ■ dU = J/" 1 -— da k 
oak 



(14) 
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of left invariant one- forms, then wedge the 15 linearly independent forms together. 1 But due to the 15 independent 
parameters needed for S'C/(4), this method is unfortunately quite time consuming and thus prohibitive. An easier 
way, initially given in ]l7p , is to calculate the 4 by 4 matrices, dll/dak (for k = {1 . .. 15}), and take the determinant 
of the coefficient matrix generated by their subsequent expansion in terms of the Gell-Mann basis. 
To begin, we take the transpose of equation (||) to generate 

u — JjT _ e «Af 5 ai5 e iAjai4 e iA3 Q13 giA^Qi2 giA^an g iA^ qio giA^ag g iA Jas 

xe ^J &7 e i^To a e £ i^3 "B e »Aj'a4 e iA3 aagiA^as e « A 3 ai Qrj\ 

An observation of the components of our Lie algebra sub-set (A2, A3, A5, As, A10, A15) shows that the transpose operation 
is equivalent to making the following substitutions 

Af — > — A 2 , — > A3, 
~ * i Ag — > Ag, 

A 10 " A 10 ) ^15 -* A 15 ( 16 ) 



in equation (US) generating 



e iAi5ai5 e iAsai4gjA3ai3 e — iA 2 ai2 e iA3Qn e — iA 5 Q!io e iA3Qci e — iA 2 Q 8 
O i\ 3 ot 7 —i\ 10 a e i\ 3 a 5 —1X5014, i\ 3 a 3 —i\ 2 a 2 „iX 3 a 



X g lA3Q7 g — lAloQ6 g lA 3 Q5 g — lA 5a4g4A 3 Q3 e — lA 2 Q2 g«A 3 Qi (17) 

Whichever form is used though, we then take the partial derivative of u with respect to each of the 15 parameters. 
In general, the differentiation will have the form 

— g« A r 5 Q i5g iA 8 "i4giA^ a 13e i\%a 12 . . , giA^a fc+1 i\ T e iX I a k g« A p . . . giAj Q l 

da k 

= g iA r5 Q i5g iA r Q i4gi A 3 r ai3gi A 2 ""i2 . . . ^"Hij^g-'^OHi . . . e ~ iX I a i2 e ~iX 3 a 13 e -i\J ai4. e -i\i 5 a 15 u ^g-j 

which, if we make the following definitions, 

C{a k ) G i * {Aj, A|\ Aj, Aj, Af , Af 5 }, (19) 

and 

£}L — gC(ai5)ai5 . . . e C(a k+1 )a k+1 
JTJ-L _ e -C(a k + 1 )a k + 1 _ _ _ g-C(ai 5 )ai 5 ^fj) 

can be expressed, in a "shorthanded" notation as 

P- = E L C{a k )E- L u. (21) 
By using these equations and the Baker-Campbell-Hausdorff relation, 

e x Ye- x = Y + [X, Y] + i [X, [X, Y}] + . . . , (22) 



we are able to consecutively solve equation (|18|) for k = {15, . . . , 1}, giving us a set of 4 by 4 matrices which can 
be expanded in terms of the 15 transposed elements of the SU(4) Lie algebra with expansion coefficients given by 
trigonometric functions of the group parameters af. 

Mk-^u- 1 = E L C(a k )E- L = J2 C ^ X J- ( 23 ) 

k 15>J>1 



1 Similarly one can wedge together the 15 right invariant one-forms which also yields the Haar measure in question. This is due to the 
fact that a compact simply-connected real Lie group has a bi-invariant measure, unique up to a constant factor. Such a group is usually 
referred to as 'unimodular' MM . 
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At this stage, we should illustrate the connection between the M^'s and the 15 left invariant one- forms that we 
could have used. To begin we note that 



du ■ u- 1 = d(U T ) ■ {U T )~ l = (dUf ■ (U-y = (U- 1 ■ dUf 
Thus the following relationship between equation (113) and (pjl) holds: 



(£ da T ^ ^ c kj Xjda k = (U-^ da k ) T . 

/ 15>j>l 



Therefore we can conclude 



dU 

C/ _1 - — dak = ^ CkjXjdak, 

"'' 1<?'<15 



(24) 



(25) 



(26) 



for k = {1, . . . , 15}. So even though we are calculating the right invariant one- forms for u, we are really calculating 
the left invariant one-forms for U. The important thing to note is that the Ckj's do not change. 2 

Now, the expansion coefficients Ckj are the elements of the determinant in question. They are found in the following 
manner 



(27) 



where the trace is done over all 15 transposed Gell-Mann matrices Jig ]. The index k corresponds to the specific a 
parameter and the j corresponds to the specific element of the algebra. Both the k and j indices run from 15 to 1. 
The determinant of this 15 by 15 matrix yields the differential volume element, also known as the Haar measure for 
the group, dVsu(4) that, when integrated over the correct values for the ranges of the parameters and multiplied by 
a derivable normalization constant, yields the volume for the group. 

The full 15 by 15 determinant Det[c/y], fc, j £ {15, . . . , 1}, can be done, or one can notice that the determinant can 
be written as 



(28) 



which differs only by an overall sign from Det [ckj] above, but which also yields a quasi-block form that generates 





Cl5,14 


Cl5,13 • 


■ Ci 5j i 


c 15,15 


Csu{i) — 


Cl4,14 


Cl4,13 • 


■ Ci 4j i 


Cl4,15 




Cl,14 


Cl,13 • 


■ c M 


Cl,15 



a 



SU(4) 






D 


A 


D 



(29) 

B is a complicated 



where D corresponds to the 9 by 9 matrix whose determinant is equivalent to dVsu(3) ■ dot\5 
6 by 9 matrix, and O is a 9 by 6 matrix whose elements are all zero. 

Now the interchange of two columns of a N by N matrix yields a change in sign of the corresponding determinant, 
but by moving six columns at once the sign of the determinant does not change and one may therefore generate a 
new matrix 



SU(A) 



Cl5,8 Cl5,7 ■ ■ • Cis,! €15,15 €15,14 €15,13 . . . €15,9 
Cl4,8 Cl4,7 ■ ■ • C144 Ci4,i5 €14,14 €14,13 ■ ■ • c 14,9 

Cj.,7 ■ ■ ■ Ci,i Ci,i5 Ci,i4 C143 . . . Ci,g 



c l,i 



(30) 



2 The transpose operation on the Gell-Mann matrices only gives an overall sign difference to some of the matrices 

f\T \T \T \T \T \T \T \T \T \T \T \T \T \T \T -1 
t A l i 1 3 ' 4 1 5 ' 6 ' 7 1 8 > 9 s A 10> A lli A 12i A 13 ' A 14: A 15/ ~ * 

{Ai, — A2, A3, A4, — A 5 , Ae, — A7, Ag, A9, — Am, An, — A12, A13, — A14, A i5 } 

but these sign changes are augmented by the inversion in the sum over k and therefore cancel out in the final construction of the 
determinant that we need. The transpose operation on equation (H) is done only to simplify the initial evaluation of the expansion 
coefficients c^j. 
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which is now block diagonal 



C'su(4) 



D 





B 


A 



(31) 



and which yields the same determinant as Cguuy Thus, with this new form, the full determinant is just equal to the 
determinant of the diagonal blocks, one of which is already known. 3 So only the determinant of the 6 by 6 sub-matrix 
A which is equal to 



C6.14 c 6,13 

C5.14 C5.13 



C6,10 C6,9 
C5,10 C5,9 



(32) 



Cl,14 Ci ; i3 . . . C140 Cl,9 

is needed. Therefore the differential volume element for SU (4) is nothing more than 

dV SU (4) = Det[cfcj] 

= -Dct[A] * Det [D]da 15 ...da>i 
= —Dct[A] * dVsu(3)dai^dae . . . dai 

which when calculated yields the Haar measure 

dVsu(4) — cos(a4) 3 cos(ae) cos(aio) sin(2a!2) sin(a4) sin(ae) 5 sin(2of8) sin(aio) 3 sin(2ai2)ciai 



. da\. 



(33) 
(34) 



This is determined up to normalization (explained in detail in Appendix |b|). Integration over the 15 parameter space 
gives the group volume 



J ■■■ J dV su(i) = (192) !■■ JdV. 



SU(4) 



SU(4) 



(35) 



which is in agreement with the volume obtained by Marinov pi 



IV. TWO QUBIT DENSITY MATRIX PARAMETRIZATION 



Using this Euler angle parametrization, any two qubit density matrix can now be represented by following the 
convention derived in Boya et. al. [js). As stated in Boya et. al., any A^-dimensional pure state can be written as a 
diagonal matrix with one element equal to 1 and the rest zero. Different classes of pure states have a different ordering 
of the zero and non-zero diagonal elements. Therefore, if one wants to write down a mixture of these different pure 
states, one must take the following convex sum 

Pd = ^ tfpu (36) 

i 

where is now the mixed state, pi (i running from 1 to N) are the pure state matrices satisfying Tr\fiiPj] = 25y, 
and the a 1 are constants that satisfy J2i a * = 1 ancl < a* < 1 @]. Now the d 1 are just the eigenvalues of the density 
matrix pd and can thus be parameterized by the squared components within the (N— 1) sphere, S 1 ^^ 1 . If we now want 
the most general mixed state density matrix in some arbitrary basis, one only has to perform a unitary, unimodular 
transformation upon pd] a transformation that will be an element of SU(N). So for our two qubit density matrix p 
we write 

P = U Pd U\ (37) 



3 The general determinant formula for this type of block matrix is given, without proof, by Robert Tucci in quant-ph\0103040. 
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where pd is the diagonalized density matrix which corresponds to the eigenvalues of the 3-sphere, S 3 ^ 0, |To[ 



Pd 



/sin 2 ( 



V 



?i)sin 2 (0 2 )sin 2 (0 3 ) 

cos 2 (6» 1 )sin 2 (6» 2 )sin 2 (6'3) 

cos 2 








9 2 )sin 2 ( 








(38) 



cos 2 (0 3 )/ 



and U, now an element of SU(A), is from equation (||). 

It is instructive to rewrite equation (|3|) as the exponentiated product of generators of the Cartan subalgebra 
that we are using in our parametrization of S77(4); e A3 * a ,e As * b , and e Al5 * c . Unfortunately, indeterminacies with the 
logarithm of the elements of pd does not allow for such a rewrite so pd will be expressed in terms of the following sum 



l<j<15 



(39) 



We begin by redefining pd in the following way 



Pd 



(w 2 x 2 y 2 

(1 - w 2 )x 2 y 2 

(1 - x 2 )y 2 

V 




(40) 



where w 2 = sin 2 (6q), x 2 — sin 2 (6> 2 ), and y 2 = sin 2 ((9 3 ). Now we calculate the decomposition of equation ( |4C| ) in terms 
of the elements of the full Lie algebra. This is accomplished by taking the trace of \pd ■ Xj over all 15 Gell-Mann 
matrices. Evaluation of these 15 trace operations yields the following decomposition of equation (MC 



p d = h i + ^(-l + 2w 2 )x 2 y 2 *\ 3 + ^(-2 + 3x 2 )y 2 *\ 8 + 1 ^=(-3 + 4y 2 )*\ 15 , 



(41) 



where the one-quarter I4 keeps the trace of pd in this form still unity. 

With equations (B^) and ( pl|) we can write down p completely in terms of the Lie algebra sub-set of the parametriza- 
tion. First, U\ the transpose of the conjugate of equation (||), is expressed as 



Jjt _ e -i^l&ais e ~i^sai4 e ~i^3a 13 e -i\ 2 a 12 g-iXzet!! e ~i\ 5 a 10 e -i\ 3 a a e -i\ 2 a8 

X g-^3<l7g-'^10«6 e -^3«5g-'^5a4g-iA3a3g-iA2O2 e -iA3ai 



(42) 



Thus equation (am is equal to 



e iA3aigiA2a2giA3Q;3 e iA5a4 e iA3a5giAioa6giA3a7 e iA2a8 

xe iA3a9 e iABQl0 e iA3ail e iA2Ql2 e iA3ai3 e iAsQl4 e iAl5ai5 

x (\U + + 2w 2 )x 2 y 2 * A 3 + ^(-2 + 3x 2 )y 2 * A 8 + ^=(-3 + 4y 2 ) * A 15 ) 

X g-iAi5ai5 e -iA 8 Qi4 e — iA 3 ai3g — iA 2 Q!i2g — iA3aiig— iAsQiOg— iA3Q9 e -iA 2 a8 
X g-iA 3 a7g — iAi Qi3g — iA3a5g-iA 5 a4g-iA3Q3g-iA 2 a2g— iA 3 ai 



(43) 



which, because I4, A3, As, and A15 all commute with each other, has the following simplification 

P = 



giA3ai3giAsai4giAi5ai5 



x (~m + + 2w 2 )x 2 y 2 * A 3 + ^=(-2 + 3s V * A 8 + -^=(-3 + Ay 2 ) * A 15 ) 

x g — iAl5Ql5 g — iA8Q!14g — iA3Q13 

• • • (il4 + + 2u, 2 ) : eV * A 3 + ^=(-2 + 3*V * A 8 + ^(-3 + 4y 2 ) * A 15 ) • 



(44) 
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Therefore, all density matrices in SU (4) have the following form 4 

p _ giA 3 ai e iA 2 a2gi- > i3a3 e iA5Q4 e iA3Q5 e iAioa 6e iA3a7 e iA2as e «^3a9 e i- > i5aiOgiA3Qii e iA2Qi2 

x (~l4 + ^(-1 + 2 W Vy 2 * A 3 + -^=(-2 + 3z V * A 8 + -^=(-3 + % 2 ) * A 15 ) 

xe -iA 2 ai2g — «A3aiig — iAsQiOg — iA3a 9 g — iA2a 8e -iA3a7 e -iAioQ;6g--iA3Q!5 e — iA 5 a4 e -iA 3 Q3 e -iA2a2 e -iA3ai (45) 

where 

w 2 = sin 2 (6»i), 
x 2 = sin 2 (0 2 ), 

y 2 = sin 2 (0 3 ), (46) 
with the ranges for the 12 a parameters and the three 9 parameters given by 

< a\, a 3 , a 5 , a 7 , ag,a n < n, 

TT 

< a 2 , a4, a e , a s , a w , a\ 2 < — , 

4 - 2' V V3 - 2'3~ 2 V ' 

where the a's are determined from calculating the volume for the group (explained in Appendix |^) and the 9's from 
generalizing the work contained in || . 

In this manner all two-particle bipartite systems can be described by a p that is parameterized using 12 Euler angles, 
and three spatial rotations, and which by (221, majorizes all other density matrices of SU(4). b Exploitations of this 
property, related to Birkhoff's theorem concerning doubly stochastic matrices and convex sets p0| , allows us to use 
this parametrization to find the subset of ranges that generate entangled density matrices and thus parameterize the 
convex polygon that describes the set of entangled two qubit systems in terms of Euler angles and spatial rotations. 
In order to do this, we need to look at the partial transpose of equation (Ha). 6 



V. REFORMULATED PARTIAL TRANSPOSE CONDITION 



To begin, one could say that a particular operation provides some entanglement if the following condition holds. 
Let p be a density matrix composed of two pure separable qubit states. Then the following matrix will represent the 
two qubit subsystems A and B, 

P = PA®PB- (48) 

Let U £ 5£/(4) be a matrix transformation on two qubits. Therefore, if 

p' = UpU^ (49) 

is an entangled state, then the operation is capable of producing entanglement p3| , B4J . One way in which we can 
tell that the matrix p' is entangled, is to take the partial transpose of the matrix and see if it is positive (this is the 
Peres-Horodecki criterion (2^, |2(|). In other words we wish to see if 

{p') TA < 0, or {p') TB < 0. (50) 

These relations imply each of the partial transposes, Ta and Tb, leaves p non-negative. If either of these conditions 
are met, then there is entanglement. 



4 One should be able to write only 12 matrices in the parametrization of U since the little group of p d is generated by A3, Ag, and A15. 

5 The eigenvalues of the given p always satisfy v\ > (1/2, V3, va) with additional ordering between the V2, V3 and eigenvalues. Therefore 
one can always find an ordering of the vi that satisfies the majorization condition. 

6 It is worth noting that Englert and Metwally fell, E2l have shown that for certain purposes, 9 parameters extracted from the density 
matrix are enough to describe certain importantcnaracteristics of the local and nonlocal properties of the density matrix. In such cases, 
one should use the density matrix representation discussed in section IVJ which more clearly expresses their ideas. 
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As an example of this we look at the situation where pd — p and U is given by equation (g) . By taking the partial 
transpose of p 1 and finding the subset of the given ranges of pd and U such that p' satisfies the above conditions for 
entanglement we will be able to derive the set of all matrices which describes the entanglement of two qubits. To do 
this, we look at the eigenvalues of the partial transpose of p' . 

Using the Euler angle parametrization previously given, a numerical calculation of the eigenvalues of the partial 
transpose of p' has been attempted. Under the standard Peres-Horodecki criterion, if any of the eigenvalues of the 
partial transpose of p' are negative, then we have an entangled p' otherwise the state p' is separable. As we have 
mentioned, we would like to derive a subset of the ranges of the Euler angle parameters involved that would yield such 
a situation, thus dividing the 15 parameter space into entangled and separable subsets. Unfortunately, due to the 
complicated nature of the parametrization, both numerical and symbolic calculations of the eigenvalues of the partial 
transpose of p' have become computationally intractable using standard mathematical software. Therefore, only a 
limited number of searches over the 15 parameter space of those parameter values that satisfy the Peres-Horodecki 
criterion have been attempted. These initial calculations, though, have shown that all possible combinations of the 
minimum and maximum values for the 12 a and three 9 parameters do not yield entangled density matrices. Numerical 
work has also shown that with this parametrization, one, and only one, eigenvalue will be negative when the values of 
the parameters give entangled density matrices. This is a verification of Sanpera et. al. and Verstraete et. al. who have 
shown that the partial transpose of an entangled two-qubit state is always of full rank and has at most one negative 
eigenvalue p7[ p8[ . This result is important, for it allows us to move away from using the standard Peres-Horodecki 
criterion and substitute it with an expression that only depends on the sign of a determinant. 

To begin, the eigenvalue equation for a 4 by 4 matrix is of the form 

(X-p 1 )(X-p 2 )(X-p 3 )(X-p i ), (51) 

which generates 

A 4 + aX 3 + bX 2 + cX + d = 0, (52) 

where 

a = -(^1+^2+^3 + ^4), 

b = fll{l2 + fJ-lfJ-3 + MlM4 + M2M3 + M2M4 + M3M4) 
C = -(/ilM2(M3 + I-M) + (/"I + P2)P3P±), 

d = (53) 

Now, since the jj,i are eigenvalues, there sum must equal 1. Thus, coefficient a in equation (|52|) is -1. Therefore the 
characteristic equation we must solve is given by 

A 4 - A 3 + bX 2 + cX + d = (54) 

which can be simplified by making the substitution r = A — 1/4 which yields 

t 4 + pr 2 + qr + r 

P 

( I 

r 

The behavior of the solutions of this equation depends on the cubic resolvent 

7 3 + 2p7 2 + (p 2 - 4r) 7 - q 2 = 0, (56) 
which has 717273 = q 2 [^9|. Recalling that the solution of a cubic equation can be obtained by using Cardano's 



0. 

b- 
b 

2 ~ 
b 

16 



3 
256' 



(55) 
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formula p9] we can immediately write down the roots for equation 

71 



72 



73 



-2p 


2i(-p 2 - I2r) 




3 


o(^zp u + z / ^ - * zj>r -+- \/ 4^— p — izr)^ -+- -\- Z(q — 


r+ 

/ zpr ) j 3 


(2p 3 


+ 27g 2 - 72pr + ^(-p 2 - 12r) 3 + (2p 3 + 27q 2 - 72pr) 2 ) i 




-2p 


323 

(1 + i\/3)(-p 2 - 12r) 




3 


Qoi ^0^3 _L 07^,2 70^^ _L / A ( <r2 1 0^3 _L ^0^3 _L 07^2 

oz 3 (^p + z iq — 1 zpr + \/ 4( — p — izr ) u + {^P + * > <7 


— i zpr ) ) 3 


(1- 


iV3)(2p 3 + 27<? 2 - 72pr + ^4(-p 2 - 12r) 3 + (2p 3 + 27g 2 - 




-2p 


62a 

(l-zV3)(-p 2 -12r) 


* 


3 


32§ (2p 3 + 27g 2 - 72pr + ^/A(-p 2 - 12r) 3 + (2p 3 + 27<? 2 


- 72pr) 2 )a 


(1 + 


iV3){2p 3 + 21 q 2 - 72pr + ^A{-p 2 - 12r) 3 + (2p 3 + 27g 2 - 


- 72pr) 2 )i 



623 

In terms of our original parameters 6, c, and d, we have for equation (t36| 



(57) 



7 3 + 2(-| - 6) 7 2 + - & + b 2 + c - 4d) 7 - 1(1 - 46 + 8c) 2 , (58) 



and therefore for its roots 



71 



= - 12(-1 + 86)(-3 + 166(1 + 6) - 16c + 64rf) + [-54(1 - 86) 2 (1 - 46 + 8c) 2 + 6^3 



x a/27(1 - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ] § 
x (12(-1 + 86) [-54(1 - 86) 2 (1 - 46 + 8c) 2 + 6^3 

x ^27(1 - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ] 3 

72 = 6(l + iV3)(-l + 86)(-3 + 166(1 + 6) - 16c + 64rf) + (-3)§ [-18(1 - 86) 2 (1 - 46 + 8c) 2 + 2\/3 

x ^27(1 - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ] i 
x (12(-1 + 86) [-54(1 - 86) 2 (1 - 46 + 8c) 2 + 6^3 

x v/27(l - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ]s)-\ 

73 = 12(1 - iy/S)(-l + 86)(-3 + 166(1 + 6) - 16c + 64d) - 3s (-3i + V3)[ - 18(1 - 86) 2 (1 - 46 + 8c) 2 + 2\/3 

x ^27(1 - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ] * 
x (24(-l + 86) [ - 54(1 - 86) 2 (1 - 46 + 8c) 2 + 6^3 

x ^27(1 - 86) 4 (1 - 46 + 8c) 4 + 16(-1 + 86) 3 (-3 + 166(1 + 6) - 16c + 64d) 3 ]3) _1 . (59) 
Now, if all three 7 solutions are real and positive, the quartic equation d55|) has the following solutions 



/7i + V72 + V73 V7i - V72 - V73 
Tl = 2 ' T2 = 2 • 



-V7i + V72 - V73 -V71 - V72 + V73 , cr ., 
T 3 = — 7, , n = 7, — ■ (60) 



Substitution of the 7 values given in equation ( |59| ) into equation ( p0| ) creates the four eigenvalue equations that the 
standard Peres-Horodecki criterion would force us to evaluate. These are quite difficult and time consuming, especially 
when 6, c, and d are written in terms of the twelve a and three 9 parameters, and can become computationally 
intractable even for modern mathematical software. But, from the previous discussion, it is obvious that with only 
one eigenvalue that changes sign, the only parameter that needs to be analyzed is d. Therefore, instead of looking at 
solutions of (pOl) one may instead look at when d from equation (M) changes sign. 7 



7 A. Wang has proposed a general solution the eigenvalue problem for the partial transpose of two qubits (see for example [ po| eq. (22)) 
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Now, the d parameter is the zeroth order A coefficient from the following equation 

Dct(p pt - I4 * A) (61) 

where p pt is the partial transpose of equation (fl5|). This is just the standard characteristic equation that yields the 
fourth-order polynomial from which the eigenvalues of the partial transpose of equation (|4^) are to be evaluated, and 
which equations (|5^) and (|5^) are generated from. Computationally, from the standpoint of our parametrization, it 
is easier to take this determinant than it is to explicitly solve for the roots of a fourth-order polynomial (as we have 
given above). The solution of equation ( pi] ) yields an expression for d in terms of the 12 a and three 9 parameters 
that can be numerically evaluated by standard mathematical software packages with much greater efficiency than the 
full Pcres-Horodccki criterion. 8 



VI. CONCLUSIONS 



The aim of this paper has been to show an explicit Euler angle parametrization for the Hilbert space of all two qubit 
density matrices. As we have stated, such a parametrization should be very useful for many calculations, especially 
numerical, concerning entanglement. This parametrization also allows for an in-depth analysis of the convex sets, 
sub-sets, and overall set boundaries of separable and entangled two qubit systems without having to make any initial 
restrictions as to the type of parametrization and density matrix in question. We have also been able to use this 
parametrization as an independent verification to Marinov's SU (4) volume calculation. The role of the parametrization 
in simplifying the Peres-Horodecki criteria for two qubit systems has also be indicated. 

Although one may generate or use other parametrizations of SU(4) and two qubit density matrices (see for example 
IpH \i% |33| , HU) our parametrization does have the advantage of not naively overcounting the group, as well as 
generating an the easily integrable Haar measure and having a form suited for generalization. Such a parametrization 
should also assist in providing a Bures distance for the space of two-qubits. Also, although previous work has been 
done on evaluatingthe eigenvalues of the partial transpose of the two qubit density matrix (for example the work 
done by Wang in p3[), our representation allows the user to effect both a reduction in the number of equations to be 
analyzed for entanglement onset from 4 to 1 while still retaining the ability to analyze the little group and orbit space 
of the density matrix as well (see for example the work contained in [^J). We also believe that this research yields the 
following possibilities: 

1. The partial transpose condition could be used to find the set of separable and entangled states by finding the 
ranges of the angles for which the density matrix is positive semi-definite. 

2. The SU(A) parametrization enables the calculation of the distance measure between density matrices and then 
use the minimum distance to a completely separable matrix as a measure of separability. Applications to other 
measures of entanglement ]3i| are straight-forward. 

3. One could use ranges of the angles that correspond to entangled states to find the ranges of the parameters in 
the parametrization in terms of the Pauli basis states by using the following parametrization for the density 
matrix 

p = jih, + a-i^i <8> h + h ® bjTj + Ckiak <8> n). (62) 
For more on this parametrization see |2l], |22j and references within. 



in which he s tates that only one equation need be evaluated to determine entanglement. Unfortunately, in order to evaluate that one 
equation (PCI] eq. (22)), six other equations must first be evaluated (pfj eqs. (23~28)). In terms of the 15 parameters needed to represent 
the Hilbert space of a two qubit density matrix, it is far easier to evaluate the zeroth order A coefficient d given in equation 1541) than 
it is to evaluate seven total equations. Even if one where to substitute and simplify, achieving one equation, its representation in terms 
of the 15 parameters needed to accurate describe the most general density matrix would still be more complicated to numerically and 
symbolic evaluate than the d parameter. 
8 This greater efficiency is based on the observation that the kernel of the mathematical software package Mathematica ver. 4.0 rel. 3, 
running on an optimized 1.5 GHz Pentium 4 Linux box with 1 gigabyte of 333MHz DDR, was unable to express equation ( pc| ) in terms 
of the 12 a and three 8 parameters in a format suitable for encoding into a C++ program. On the other hand, it was quite easy to 
obtain all the coefficients of equation ( [54] ) in terms of our Euler parameters, simplify them, and encoded them into a CH — h program for 
numerical evaluation. Also, since only one eigenvalue ever goes negative, we have reduced the number of equations to solve from 4 to 1 
which is a definite improvement in calculatory efficiency. 
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4. Related to this last question is the question of the boundary between the convex set of entangled and separable 
states of the density matrices. For example one could use the explicit paramctrization to calculate specific 
measures of entanglement like the entanglement of formation for different density matrices in different regions 
of the set of density matrices and see which regions of the convex set correspond to the greatest entanglement 
of formation. Another possibility is that given the boundary in the a, t form, we could recreate it in terms of 
the Euler angles. 

There are obviously more, but for now, it is these areas that we believe offer the most interest to those wishing to 
develop a deeper understanding of bipartite entanglement. Also, since the methods here are quite general and rely 
primarily on the group theoretical techniques developed here, we anticipate generalizations to higher dimensional 
state spaces will be, in principle, straight-forward. 
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APPENDIX A: COMMUTATION RELATIONS FOR SU(4) 
We first note that the Gell-Mann type basis for the Lie algebra of SU (4) is given by the following set of matrices 



Ai = 



A 4 = 



A 7 = 



Aio — 



U3 




A 2 = 



A s = 



i 


Vo 



-i 






As = 



i 



/ -i 



1 

Vo 

fl 
1 


Vo o 

/° 



An = 









-2 






1 


Vo i o o 

/0 



-i 

Vo i 



A 3 = 



A fi = 



Al2 — 



A IF 




(Al) 



In order to develop the Cartan decomposition of SU(4) it is helpful to look at the commutator relationships between 
the 15 elements of its Lie algebra. In the following tables we list the commutator solutions of the corresponding zth 
row and jth column Gell-Mann matrices corresponding to the following definitions 

[Aj, Xj] =2ifijkXki 

fijk =^- i Tr i[K,X j } ■ A fe ]. 



Table 1 : [kufa] G L(K) 

(on the following page) This table corresponds to the L(K) subset of SU(4), {Ai, . . . , Ag, A15} and shows that 
for kuh &L{K), [fa,k 2 ] eL(K). 

Table 2 : [pi,p 2 ] £ L{K) 

(on the second page) This table corresponds to the L(P) subset of SU(A), {Ag,...,Ai4} and shows that for 
Pi,P2 € L(P), [ Pl ,p 2 ] G L(K). 

Table 3 : [ki,p 2 ] G L(P) 

(on the third page) This table corresponds to the commutator solutions for the situation when k\ G L(K) 
andp 2 G L(P), [ki,p 2 ] G L(P). 



TABLE I: [ki,k 2 ] € L(K) 





Ai 


A 2 


A 3 


A 4 


A 5 


A 6 


A7 


As 


A15 


Ai 





2iA 3 


-2iA 2 


iA7 


— iA6 


iAs 


— iA 4 








A 2 


-2iA 3 





2«Ai 


iA 6 


zA7 


— iA 4 


— 1A5 








A 3 


2iA 2 


-2iAi 





iA 5 


— iA 4 


— jA7 


iA 6 








A 4 


— i\y 


— iX e 


— iA 5 





*(A 3 + V^3A S ) 


iA 2 


iAi 


-iV^As 





A 5 


i\e 




iA 4 


-i(A 3 + V^3A 8 ) 





— iAi 


iA 2 


iv^A^i 





Ae 


-iA 5 


iA 4 


iA7 


— iA 2 


iAi 





»(-A 3 + v^Ag) 


— i V3A7 





At 


iA 4 


iAs 


— iAe 


—iAi 


— iA 2 


j(A 3 - V3A 8 ) 





iv^Ae 





A 8 











iV3 A 5 


— i V3A 4 


i\/3A7 


— i\/3A6 








Al5 



























TABLE II: \pi,p 2 ] € L(K) 





Ag 


Am 


An 


A12 


A13 


A14 


A 9 





i(A 3 + ^A s + 2y§A 15 ) 


iA 2 


2A1 


iA 5 


iA4 


Am 


-i(A 3 + ^- A 8 +2^/|Ai 5 ) 





-»Ai 


i\ 2 




iAs 


An 


-1A2 


»Ai 





J (-A3 + ^=A 8 + 2^/|Ai 5 ) 




iA 6 


A12 


— zAi 


— iA2 


i(A 3 -^A 8 -2^/|Ai5) 





— iAe 


2A7 


A 13 


— iA 5 


iA4 


— i\f 


iA 6 







Al4 


— iA4 


— iA 5 


-iXe 


— 1X7 








TABLE III: [fci,p 2 ] € L(P) 





Ai 


A 2 


A 3 


A 4 


As 


Ae 


At 


As 


A 9 


Aio 


An 


Al2 


Al3 


Al4 


A 15 


Ai 


















»Ai2 


— iAn 


iAio 


— iXg 










A 2 


















iAn 


iAi2 


-*A 9 


— iAio 










A 3 


















iAio 


— iXg 


-»Al2 


iAn 










A 4 


















iAi4 


— »Al3 








iAio 


-iAg 




As 


















iAl3 


iAi4 








— iAg 


-iAio 




Ae 
























lAl4 


-iAi 3 


iAi2 


-iAn 




At 
























»Al 3 


iAi4 


— iAn 


-iAi2 




As 


















7f A1 ° 


"73 Ag 


73 Al2 


"73 Al1 




73 A - 




A 9 


-iAi2 


— »An 


—iAio 


— »Al4 


-iAi 3 






















—«^/ f Aio 


Aio 


iAn 


— »Al2 


iAg 


iAi 3 


— ?Ai4 








7a A ° 














i\/§A 9 


An 


—iAio 


iAg 


iAi2 








— »Al4 


-»Al3 


~vI Al2 
















A12 


iAg 


iAio 


— »An 








iAi3 


— iAi4 


75 Al1 














V I A n 


A13 











— iAio 


iAg 


-«Ai2 


iAn 


75 Al4 














— i-^ §Ai4 


Al4 











iXg 


iAio 


iAn 


iAi2 


~73 Al3 














V 3 Al3 


A15 




















~V§ A9 


V§Ai2 


"V |An 


V |Al4 


-*y§Ai 3 
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APPENDIX B: INVARIANT VOLUME ELEMENT NORMALIZATION CALCULATIONS 



Before integrating dVsuU) we need some group theory. We begin with a digression concerning the center of a group 
p6| . [37j. If S is a subset of a group G, then the centralizer, Cq(S) of S in G is defined by 

C(S) = C G (S) = {x £ G I if s € S then xs = sx}. (Bl) 

For example, if S — {y}, C(y) will be used instead of C({y}). Next, the centralizer of G in G is called the center of 
G and is denoted by Z(G) or Z . 

Z(G) = Z = {z e G I zx = xz for all x £ G} 

= c o(G ). <B2) 

Another way of writing this is 

Z(G) = n{C(x) \xeG} 

= {z I if x £ G then z £ C(x)}. 1 ' 

In other words, the center is the set of all elements z that commutes with all other elements in the group. Finally, 
the commutator [x, y] of two elements x and y of a group G is given by the equation 

[x,y] = x~ 1 y~ 1 xy. (B4) 

Now what we want to find is the number of elements in the center of SU(N) for N = 2, 3, and 4. Begin by defining 
the following 

Z n = cyclic group of order n = Z„ = Z(SU(N)). (B5) 

Therefore, the set of all matrices which comprise the center of SU(N), Z(SU(N)), is congruent to Zjv since we know 
that if G is a finite linear group over a field F, then the set of matrices of the form T*c g g, where g £ G and c g £ F, 
forms an algebra (in fact, a ring) [|l5| [37|. For example, for SU(2) we would have 

Z 2 ={x £ SU{2) I [x,y] £ Zi for all y £ SU(2)}, 
[x,y]=ul a , (B6) 
Zi ={1 2 }. 

This would be the set of all 2 by 2 matrix elements such that the commutator relationship would yield the identity 
matrix multiplied by some non-zero coefficient. In general this can be written as 

Z N ={x £ SU(N) I [x, y] £ Zi for all y £ SU(N)}, 

Z X ={*n}. (B7) 

This is similar to the result from |36|, that shows that the center of the general linear group of real matrices, GL^(JR), 
is the group of scalar matrices, that is, those of the from col, where I is the identity element of the group and u> is 
some multiplicative constant. For SU(N), uil is an N th root of unity. 

To begin our actual search for the normalization constant for our invariant volume element, we first again look at 
the group 577(2). For this group, every element can be written as 

(B8) 




where \a\ 2 + \b\ 2 = 1. Again, following [^6| we can make the following parametrization 

a = yi- iy 2 , 

b = y 3 -iyi, (B9) 

1 = 1/1+02 + 1/3 + vl 

The elements (1,0,0,0) and (—1,0,0,0) are anti-podal points, or polar points if one pictures the group as a three- 
dimensional unit sphere in a 4-dimensional space parameterized by y, and thus comprise the elements for the center 
group of SU(2) (i.e. ±12). Therefore, the center for SU{2) is comprised of two elements. 
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In our parametrization, the general SU(2) elements are given by 

dVsu(2) — sia(2v) dfidud^, 



with corresponding ranges 



0</i,£<7T, 
7T 

0<v<~. 



(BIO) 

(BH) 
(B12) 



Integrating over the volume element dVsu(2) with the above ranges yields the volume of the group SU(2)/ 'Z 2 - In other 
words, the SU(2) group with its two center elements identified. In order to get the full volume of the SU(2) group, 
all ones need to do is multiply the volume of SU(2)/Z2 by the number of removed center elements; in this case 2. 

This process can be extended to the SU(3) and SU(4) parametrizations. For SU(3) J6|, §, |9[ |l0| (here recast as a 
component of the ST/ (4) parametrization) 



SU(3) = e iA3a7 e iA2a8 e iA3a9 e lA5Ql0 D(aii,ai2,ai 3 )e lA8C 



(B13) 



Now, we get an initial factor of two from the D(qu, ai2, 013) component. We shall now proves that we get another 
factor of two from the e lA3a9 e lA5ai ° component as well. 

From the commutation relations of the elements of the Lie algebra of SU(3) (see || for details) we see that 
{A3, A4, A5, As} form a closed subalgebra SU{2) x f7(l). 9 

[A3, A 4 ] =i\ 5 , 
[A3, A 5 ] = - «A 4 , 
[A 3 ,A 8 ] =0, 

[A 4 ,A 5 ] =i(A 3 + V3A 8 ), 
[A 4 , A 8 ] = - iV3X 5 , 

[A 5 , A 8 ] =iV3X4. (B14) 

Observation of the four A matrices with respect to the Pauli spin matrices of SU(2) shows that A 4 is the SU(3) 
analogue of cti, A5 is the SU(3) analogue of cr 2 and both A3 and A 8 are the SU(3) analogues of (73 



01 



(T-2 



0-3 




A 4 = 




(B15) 



and A 8 = 




\0 -2j 



Thus one may use either {A3, A5} or {A3, A5, A 8 } to generate an SU{2) subgroup of SU{3). The volume of this SU{2) 
subgroup of SU(3) must be equal to the volume of the general SU(2) group; 2n 2 . If we demand that any element of 
the SU(2) subgroup of SU(3) have similar ranges as its SU{2) analogue 10 , then a multiplicative factor of 2 is required 
for the e iX 3 a s e iX 5 a 10 com p 0n ent. 1:L 



9 Georgi [ p8[ has stated that A2, A5 , and A7 generate an SU (2) subalgebra of SU (3). This fact can be seen in the commutator relationships 
between these three A matrices contained in |p| or in Appendix M 

10 This requires a normalization factor of — j= on the maximal range of As that is explained by the removal of the Z3 elements of SU (3) . 

11 When calculating this volume element, it is important to remember that the closed subalgebra being used is SU(2) X (7(1) and therefore 
the integrated kernel, be it derived either from e tXaa e lA 5' 3 e lA37 or e tX3a e tX5 @ e lX8 ~ l , will require contributions from both the SU(2) and 
U(l) elements. 
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Finally, SU(3) has a Z$ whose elements have the generic form: 














m 





V° 








where 

vl = vl = 1. 

Solving for 7?i and 772 yields the following elements for Z3 




(-1) 














1 




■-( 













(-1) 








(-1)1 








(-1) 



(B16) 



(B17) 



(B18) 



which are the three cube roots of unity. Combining these 5/7 (3) center elements, a total of three, with the 2 factors 
of 2 from the previous discussion, yields a total multiplication factor of 12. The volume of SU(3) is then 

V su(3) = 2*2*3*V(SU(3)/Z 3 ) 

= V3tt 5 (B19) 

using the ranges given above for the general 5/7(2) elements, combined with < au < -^=. Explicitly: 

< a 7 , a 9 ,an,a 13 < n, 
< a 8 ,a w ,ai2 < —, 

< < -^=. (B20) 
v3 

These are modifications of || pi |l(| |3{| and take into account the updated Marinov group volume values p9| . 

For 5/7(4) the process is similar to that used for 5/7(3), but now with two 5/7(2) subgroups to worry about. For 
5/7(4), 



U 



[SU(3)]e 



iAi 5 Qi5 



(B21) 



Here, the two 5/7(2) subalgebras in 5/7(4) that we are concerned with are {A3, A4, A 5 , As, A15} and {A3, A9, Ai , Ag, A15}. 
Both of these 5/7(2) x Z7(l) x /7(1) subalgebras are represented in the parametrization of 5/7(4) as 5/7(2) subgroup 
elements, e 4A3a3 e iA5 " 4 and e ,A3 " 5 e ,Jl0 " e . We can see that A 10 is the 5/7(4) analogue of cr 2 12 and A i5 is the 5/7(4) 
analogue to 03 13 . The demand that all 5/7(2) subgroups of 5/7(4) must have a volume equal to 27r 2 is equivalent to 
having the parameters of the associated elements of the 5/7 (2) subgroup run through similar ranges as their 5/7 (2) 
analogues. 14 As with 5/7(3), this restriction yields an overall multiplicative factor of 4 from these two elements. 15 
Recalling that the 5/7(3) element yields a multiplicative factor of 12, all that remains is to determine the multiplicative 
factor equivalent to the identification of the 5/7(4) center, Z4. 

The elements of the center of 5/7(4) are similar in form to the ones from 5/7(3); 



(B22) 



^1 








\ 





172 














'/3 





\o 











We have already discussed A5 in the previous section on SU(3). 
It is the SU (4) Cartan subalgebra element. 

This requires a normalization factor of -4= on the maximal range of A15 that is explained by the removal of the Z4 elements of 5(7(4). 
When calculating these volume elements, it is important to remember that the closed subalgebra being used is SU(2) X U(l) X (7(1) 
and therefore, as in the SU(3) case, the integrated kernels will require contributions from appropriate Cartan subalgebra elements. 
For example, the e 1 ^ 3 " 3 e'* 5 ™ 4 component is an 5(7(2) sub-element of the parametrization of 5(7(4), but in creating its corresponding 
5(7(2) subgroup volume kernel (see the S!7(3) discussion), one must remember that it is a 5(7(2) C 5(7(3) C 5(7(4) and therefore the 
kernel only requires contributions from the A3 and \g components. On the other hand, the e lA3Q5 e lAl0 " 6 element corresponds to a 
5(7(2) C 5(7(4) and therefore, the volume kernel will require contributions from all three Cartan subalgebra elements of 5(7(4). 
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where 

Vt = Vl =4 = 1- (B23) 

Solving yields the 4 roots of unity: ±14 and ± H4, where I4 is the 4x4 identity matrix. So we can see that Z4 
gives another factor of 4, which, when combined with the factor of 4 from the two SU (2) subgroups, and the factor 
of 12 from the SU (3) elements, gives a total multiplicative factor of 192. Integration of the volume element given in 
equation ([54]) with the following ranges 

< a\, a 3 , a 5 , a 7 , ag, a n , "13 < n, 

7T 

< a%, 0*4, a e ,a$, a w , an < — , 
0<a 14 <^, 

< aiB < (B24) 

gives 

V SU {4,) = 2*2*2*2*3*4* ^(S , C/(4)/Z 4 ) 

^ (B25) 



3 

This calculated volume for SU (4) agrees with that from Marinov |l 
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APPENDIX C: MODIFIED PARAMETER RANGES FOR GROUP COVERING 

In order to be complete, we list the modifications to the ranges given in Appendix [b] that affect a covering of SU(2), 
SU(3), and SU(4) without jeopardizing the calculated group volumes. 

To begin, in our parametrization, the general SU(2) elements are given by 

dVsuti) = sin(2i/) d/j,dvd£, 
with the corresponding ranges for the volume of SU(2)/Z 2 given as 

0</U,£<7r, 

<i/<|. (C2) 

In order to generate a covering of SU(2), the £ parameter must be modified to take into account the uniqueness of 
the two central group elements, ±12, under spinor transformations. 16 This modification is straightforward enough; 
£'s range is multiplied by the number of central group elements in SU{2). The new ranges are thus 

< ft < it, 

IT 

0<u<-, 

< £ < 2tt. (C3) 

These ranges yield both a covering of SU(2), as well as the correct group volume for SU{2). 17 

For SU(3), here given as a component of the SU(4) parametrization, we know we have two SU{2) components 
(from Appendix |b|) , 

SU(3) = e iX3a7 e iX2aa e iX3a9 e iX5ai0 D(a n ,ai2, a 13 )e lXsaii . (C4) 

Therefore the ranges of ag and a\ 3 should be modified just as £'s was done in the previous discussion for SU(2). 
Remembering the discussion in Appendix |b| concerning the central group of 577(3), we can deduce that a 14 's ranges 
should be multiplied by a factor of 3. This yields the following, corrected, ranges for SU(3) 18 

< «7, an < 7r, 

7T 

< a 8 , aio, ai2 < —, 
< ag, ai3 < 27r, 

< a u < V3ir. (C5) 

These ranges yield both a covering of SU(3), as well as the correct group volume for SU(3). 
For SU(4), we have two SU{2) subgroup components 

SU(4) = e ^3ai e <A2a2 e iA3a3 e iA 5 a4 e a 3 a 5e iAioa 6 [ iSU -(-3- ) ] e iAi5ai5_ ( C6 ) 

As with the SU{2) subgroup ranges in SU(3), the ranges for a 3 and as each get multiplied by 2 and ais's ranges get 
multiplied by 4 (the number of SU (4) central group elements) . The remaining ranges are either held the same, or 



For specific examples of this, see either J12| or . 

One may interchange fi and £'s ranges without altering either the volume calculation, or the final orientation of a two-vector under 
operation by D. This interchange is beneficial when looking at Eulcr para,mct r iza t ions beyond SU(2). 

Earlier representations of these ranges for SU (3). for example in [H, H, H, lloL tq, B9| , were incorrect in that they failed to take into 
account the updated SU (N) volume formula in ]19| . 
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modified in the case of the SU(3) element; 

< ai, 07, an < t, 
< a 2 , a4, "6, "8, «io, «12 < ^, 
< a 3 ,a 5 , a 9 , a i3 < 2-k, 
< aw < V3tt, 

< aie < 2^/|tt. (C7) 

These ranges yield both a covering of SU{ / \) 1 as well as the correct group volume for SU(A). 

In general we can see that by looking at SU(N)/Zn not only can we arrive at a parametrization of SU(N) with a 
logically derivable set of ranges that gives the correct group volume, but we can also show how those ranges can be 
logically modified to cover the entire group as well without any arbitrariness in assigning values to the parameters. 
It is this work that will be the subject of a future paper. 
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